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Abstract 

We have studied the monopole-percolation phenomenon in the four dimen- 
sional Abelian theory that contains compact U(l) gauge fields coupled to uni- 
tary norm Higgs fields. We have determined the location of the percolation 
transition line in the plane (/3g,/3ff). This line overlaps the confined-Coulomb 
and the confined-Higgs phase transition lines, originated by a monopole- 
condensation mechanism, but continues away from the end-point where this 
phase transition line stops. In addition, we have determined the critical expo- 
nents of the monopole percolation transition away from the phase transition 
lines. We have performed the finite size scaling in terms of the monopole 
density instead of the coupling, because the density seems to be the natural 
parameter when dealing with percolation phenomena. 
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I. INTRODUCTION 



Much work has been devoted in recent years to the understanding of the pure gauge U(l) 
lattice field theory. Despite the apparent simplicity of this model, which can be considered 
a limit of the well known Z(n) gauge models in four dimensions, the order of their phase 
transition is a rather controversial issue. Pioneering simulations suggested that the transition 
was actually of second order, i.e., a continuum limit for this lattice theory is possible [|TJ. 
Posterior analysis showed the presence of a two-peak structure in the histograms questioning 
the order of the phase transition Several new approaches have been proposed: extended 
lattice actions ||, different topologies monopole suppressed actions [0], etc. For a 

summary of the history of this transition see Ref. ||. 

The main point in order to clarify the nature of the phase transition in pure gauge 
U(l) lattice field theory is the understanding of the role of topological excitations. This 
analysis was initiated in the early eighties by several groups: Einhorn and Savit and 
Banks, Kogut and Myerson [|Hj. Their conclusion was that topological excitations are 
actually strings of monopole current. In addition, DeGrand and Toussaint showed that 
these monopoles can be actually studied by direct numerical simulation, in particular the 
monopole condensation phenomenon over the phase transition. On the other hand, it is 
widely accepted that the observed problems in the analysis of the phase transition can be 
related to the one-dimensional character of the topological excitations. 

An important observation was made by Kogut, Kocic and Hands ]P2]] . They showed that 
in pure gauge non-compact U(l) (i.e. the action obtained keeping only the first order in 
the Taylor expansion) which is Gaussian, monopoles percolate and satisfy the hyperscaling 
relations characteristic of an actual second-order phase transition. 

Baig, Fort and Kogut |L3| showed that in the compact pure gauge theory, just over the 
phase transition point, monopoles not only condensate but also percolate. They pointed out 
that the strange behavior of this phase transition -its unexpected first order character and 
strange critical exponents- can be related to the coincidence of these two phenomena. 

Furthermore, Baig, Fort, Kogut and Kim [14| showed that in the case of non-compact 
QED coupled to scalar Higgs fields, the monopole percolation phenomena -previously ob- 
served over the gauge line- actually propagate into the full (P 9 ,Ph) plane, being (3 g and 
the gauge and Higgs couplings. 

It is surprising that the monopole percolation phenomenon is not related to the phase 
transition line that separates the confined and the Higgs phases, a second order transition 
that is logarithmically trivial. 

In the present paper we have performed a numerical simulation of the compact lattice 
gauge U(l) coupled to Higgs fields of unitary norm. First of all, we have reproduced some 
results concerning the phase diagram previously obtained by Alonso et al. |15|], but we 
have measured at the same time the behavior of monopoles, condensation and percolation. 
We have determined the location of the percolation transition line -the line defined by the 
maxima of the monopole susceptibility- in the plane (f3 g , fin)- We have observed that this line 
overlaps the confined- Coulomb and the confined- Higgs phase transition lines, which result 
by a monopole-condensation mechanism. But it continues away from the end point where 
this phase transition line stops. In addition, we have determined the critical exponents of the 
monopole percolation transition away from the phase transition lines. We have performed 
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the finite size scaling in terms of the monopole density instead of the coupling, because the 
density seems to be the natural parameter when dealing with percolation phenomena. 

The plan of the paper is the following: In Sec. [TT] we define our model and we analyze 
the appearance of the topological excitations. In Sec. [I I J we extend the standard perco- 
lation analysis to the case of monopole percolation, using densities as critical parameters. 
Sec. [IV A| contains the result of the relation between monopole percolation and phase tran- 
sition Sec. [IV B| is devoted to the finite size analysis of the percolation phenomena. Finally, 
Sec.[V| contains the conclusions of this work.[] 



II. THE MODEL 
A. The action 



In this paper, we resume the work initiated in Ref. [13], and focus on the compact abelian 



Higgs Model with a fixed length scalar field whose action is 

S = ~Pg cos W _ Ph {<f>l u x,n<f>x+n + cc.) , (1) 

where O^n) is the circulation of the compact gauge field around a plaquette, j3 g = 1/e 2 is 
the gauge coupling and 4> x = expia(x) is a phase factor. We choose fixed length Higgs fields 
because it has been shown that they live in the same universality class than conventional 
variable length Higgs fields but have the advantage that they do not require fine-tuning W7 



In addition, we want to keep as close as possible to the analysis of the non-compact Abelian 
Higgs Model of Refs. |14],18[1 in order to compare to non-compact scalar electrodynamics, 
whose action is 

Sue = \p 9 £ e, v (n) 2 -p H J2 (€U x ^ + » + cc.) . (2) 

* nfiv x,fj. 

The phase diagrams for the two models are well known. They are qualitatively represented 
in Fig. [l|. Non-compact scalar QED has two disconnected phases: a Coulomb phase and 
a Higgs phase. It is important to remark that there is not phase transition over the pure 
gauge axis. In contrast, in the compact case a transition separates the Coulomb and the 
confined phase while the transition separating the confined and Higgs phases has an end 
point (point E in the Fig. |TJa) . 



1 While our work was being completed, it appeared a work by Franzki, Kogut and Lombardo [16] 
where they analyze a model with both, scalars and fermions. In the conclusions we compare the 
results coming from both simulations and, in particular, we discuss the different ways of extracting 
the critical exponents. 
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(a) (b) 
FIG. 1. Qualitative phase diagrams for the compact and non-compact scalar QED. 

B. Monopoles 



Following |TIJ| we can separate the plaquette angle into two pieces: physical fluc- 
tuations which lie in the range — 7r to tt and Dirac strings which carry 2tt units of flux. 
Introducing an electric charge e we define an integer-valued Dirac string by 



eO 



(3) 



where the integer S^ v determines the strength of the string threading the plaquette and 
eQ^ lu represents physical fluctuations. The integer- valued monopole current, m^n), defined 
on links of the dual lattice, is then 



1 

— < 

2 



in j n ) -c^aA+S^aO + p>), 



where A 



(4) 



is the oriented sum of the 



is the forward lattice difference operator, and m M 
^ v around the faces of an elementary cube. This gauge-invariant definition implies the 
conservation law A~m At (fi) = 0, which means that monopole world lines form closed loops. 

The appearance of topological excitations in the abelian Higgs model was studied in 
Ref. |9|,|TIJ. Furthermore, in Ref |19| a Monte Carlo simulation was performed to look for 
the density of monopoles. They observed that the confined (Higgs and Colulomb) phase 
was characterized by a large (vanishing) density of monopoles vortices and electric current 
densities. 



III. PERCOLATION 

A. Concepts of percolation 

The role of the monopole percolation in this model is analyzed using the techniques of 
standard percolation EQ|. In the simplest models the sites are occupied with a probability p. 
One can define a cluster as a group of occupied sites connected by nearest-neighbor distances. 
When a cluster becomes infinite in extent, it is called a percolation cluster. Obviously, for 
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p = all the sites are empty and there is no percolation cluster. For p = 1 all the sites are 
occupied and one percolation cluster exists. There exists a critical concentration p c such 
that for p < p c (p > p c ) no (one) percolation cluster exists. This is the typical behavior of 
a phase transition. We define as the probability for an occupied site to belong to the 
infinite cluster. 

If we define g n as the relative number of clusters of size n, the probability of any site to 
belong to a cluster of size n is ng n . According to this definitions, the mean cluster size is 

s= Enn 2 g n ^ (5) 

where the infinite cluster is excluded from the sum. This quantity diverges at the critical 
concentration. 

Thus, Pqo is the order parameter of the transition and S is its associated susceptibility. 
Their behaviors near the critical point are: 

Poo ~ (p - Pc) 13 for p > Pc, (6) 

and 

s~\p- Pe r. (?) 



B. Percolation of monopoles in QED 

We have defined the monopole current m^x) in eq. (Q). One can define a connected 
cluster of monopoles as a set of sites joined by monopole line elements [^TJ. Notice that this 
definition ignores the fact that monopoles are actually vectors. The density of occupation is 

Umax 

P = T? = £ 9nn, (8) 

L n=4 

where n to t is the total number of connected sites and n max is the size of the largest cluster. 
The number 4 in the sum comes from the conservation law. 
The order parameter (Poo) is 

M = (9) 

n to t 

Its corresponding susceptibility (the mean cluster size S) is 

X = — [Yl 9 n n 2 - n 2 max . (10) 



n tot 



n=4 



Near the critical point, M should have a non analytical behavior with a "magnetic" 
exponent /3. We have followed lines of constant @h, so 

M ~(A C -A/ for ( n ) 
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where the upperscript c means "critical" . The susceptibility also diverges: 

x~(A c -<T- (12) 

One should remark that the natural parameter when studying percolation is the proba- 
bility of occupation (in random percolation it determines all the observables) . In a lattice 
gauge model, the fundamental variables are the couplings, which determine all the other 
observables. In particular, they determine the density of monopoles p, i.e., the probability 
of occupation. In this sense, it is also possible to parameterize percolation as a function of 
the densities. In this case the critical behavior will be 

M~(p-p c f for p g <p c g , (13) 

and 

X~(p-PoP. (14) 

Clearly, the two parameterizations might have the same critical exponents, but the first one 
seems to exhibit stronger finite size effects. 



IV. NUMERICAL RESULTS 

A. Phase transition vs. monopole percolation 

We have performed numerical simulations with the action (|l|) on hypercubical lattices 
with standard periodic boundary conditions. In these simulations we have measured the 
internal gauge and Higgs energies 

E g = 53 cos (^(n), (15) 

n/j,v 

E H = J2(€U x ,^ x+ , + c.c). (16) 

as well as all magnitudes related to monopole percolation, i.e. p, M and \, which are defined 
in eqs. (@, ||, [L(|) respectively. 

In order to locate the continuation into the full (/3 g , (3h) plane of the monopole percolation 
transition previously determined over the pure gauge axis in Ref |13|], we have performed 
repeated measurements over thermal cycles changing the (3 g coupling, until the entire region 
inside the limits < (3 g < 1.5 and < /3h < 1-5 has been covered. In particular, the confined- 
Coulomb and confined-Higgs phase transition lines. Lattice size for this exploratory runs 
has been of 6 4 and the j3 g step has been A/3 g = 0.005. 

We have not used any reweighting technique when measuring the monopole-related mag- 
nitudes. In all the plots, lines are only to guide the eyes and the errors are not shown because 
they are smaller than the symbols used. 

The number of iterations at each f3 g value has been 30.000, discarding the first 10.000. 
In addition, several thermal cycles have been performed at different values of j3u in order to 
determine the location of the Coulomb- Higgs phase transition line. 
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FIG. 2. The plane ((3 g , 0h) showing the lines of phase transition and the monopole percolation. 
Thermal cycles have been done at fixed 0h values: A {(3h = 0.00), B = 0.25), C ((3h = 0.30) 
and D [f3u = 0.30). The magnitudes measured along this lines are shown in the Figs, ^a, |||b, ||c 
and f4|d respectively. The point E is the end of the phase transition. 



Fig. ^ collects the results of these runs. The small black points represent the maximum 
of the monopole susceptibility. The solid line is the phase-transition line. Note that the 
monopole percolation transition line is on top of both the confined- Coulomb line and the 
confined-Higgs line up to the end point. Beyond this point monopole percolation transition 
line continues, now unrelated to any energy phase transition, approaching the vertical axis. 
This behavior is in clear contrast with the non-compact case |14| where the monopole- 
transition is completely independent of the Coulomb-Higgs phase transition. 

In Fig. |3] we compare two sets of histograms corresponding to cycles: (a) for Ph = 0.25, 
over the confined-Higgs phase transition line and (b) for Ph = 0.3, i.e. just beyond the 
end-point E. In the two figures the solid line corresponds to the value of (3 g for which the 
monopole susceptibility has a maximum. The appearance of a two-peak structure over the 
phase transition confirms that the lattice size, the statistics used and the coupling step in 
the thermal cycle are enough to establish the appearance of the phase transition. 

An interesting question to establish is the concordance of all the phenomena (phase 
transition, monopole condensation and monopole percolation) at the same coupling. In this 
sense it was proposed in Ref. |13|] to use M as an order parameter to determine the location 
of the phase transition. The existence of an end-point for the confined-Higgs phase transition 
line gives us a chance to check the behavior of these three parameters over and right beyond 
the end-point. In Fig. ^we have collected the results for thermal cycles at different values of 
Ph- Fig- corresponds to the pure gauge compact action (point A in the phase diagram), 
the case of Ref. JT3|. Fig. [|b corresponds to Ph = 0.25, i.e., over the confined-Higgs phase 
transition (point B). The behavior of all three parameters is similar to that observed in 
the pure gauge case, although the discontinuity is more abrupt. Fig. f|c corresponds to 
Ph — 0.30 (point C), just after the end point, and Fig. f|d to Ph = 0.5 (point D), a line 
that crosses the maximum of the monopole susceptibility but is far away from the end-point. 
In these last cases the energy clearly shows no discontinuity and the monopole density is 
completely smooth. Nevertheless, the M parameter shows a fast decrease corresponding to 
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6 , P H =0.25, 10.000+20.000 it. 



6 , (5 H =0.3, 10.000+20.000 it. 
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P,=0.785 




c o p s =0.790 




0.50 



(a) (b) 
FIG. 3. Some histograms. The plot (a) is evaluated in the point B. Notice the clear two peaks 
structure revealing a first order phase transition. This structure is not present in the plot (b), 
evaluated in the point C where no specific heat peak is present. 



the percolation phenomena investigated in the next subsection. 

Fig. H] shows another hint of the change of nature of the percolation at the end-point. 
We show the maximum of \ along the percolation line. A clear change can be seen at the 
end point. 



B. Finite size scaling of monopole percolation transition 

We want to determine the critical exponents of the monopole percolation away from the 
phase transition lines, at the point D = (0.5,0.606) in the plane (P 9 ,Ph)- To do this, we 
keep Ph — 0.5 fixed, and we vary (3 g within the range [0.594, 0.620] for lattice sizes ranging 
from 8 4 to 12 4 . We have recorded some percolation parameters for each simulation: the 
density of monopoles p, M and x- The statistics for each point is of 10.000 iterations once 
the first 5.000 have been discarded. 

In Fig. H we collect the results for M and x as a function of j3 g . In Fig. |7|, we show 
the same results for M and x as a function of p. 

According to finite size scaling arguments |22| , the peak of the susceptibility should grow 
with the lattice size L as 

Xmax ~ (17) 

and the value of M at the critical point must vanish as 

M crit ~ L-V". (18) 

We show the value of Xmax as a function of L in a log-log plot in the Fig. |a. Notice that 
the scaling law is satisfied. A fit to eq. (|i~7|) gives 

- = 2.38(1). (19) 

The straight line resulting from the fit is also shown. 
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Susceptibility peaks along the percolation line 
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FIG. 5. Susceptibility peaks along the percolation line. Notice the change of behavior at the 
end of the phase transition. 
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(a) (b) 
FIG. 6. M (in a) and x (i n b). Both as a function of /3 9 , keeping (5h 
effects of M are stronger than in Fig. [7|a. 



0.5. The finite size 



Scaling of M(p) (P H =0.5) 



Scaling of x(p) (P„=0.5) 
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(a) (b) 
FIG. 7. M (in a) and x (i n b). Both as a function of p, keeping (3h = 0.5. The finite size 
effects of M are lesser than in Fig. ya. 
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FIG. 8. Log-log plots for Xmax(L) and M(p — p c ), being p c = 0.377. 
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We measure (3 fitting the values that are not distorted by finite-size effects to the equa- 
tion (O). The result is 



P = 0.483(4). 



(20) 



In Fig. |b we show the selected points and the fit, with a value for the critical density about 
of p c = 0.377. 

The fit of (|TJ) is not very good. It gives j3/v = 0.7(1). Putting our measurements of 
j/u into the hyperscaling relation 



I 
v 



we obtain 







v 



0.81(1), 



(21) 



(22) 



which is more accurate than our measurement. So, we take this value and regard the 
measurement as a test of hyperscaling. 



Combining eqs. |20] and the value of v is 

v = 0.60(1). 



(23) 



The remaining critical exponents can be derived from those above using hyperscaling 
relations. 

In the below table we compare our results with other related works: 





Pc 


P 


V 


Our results 


0.377 


0.483(4) 


0.60(1) 


Non compact quenched QED |12| 




0.58(2) 


0.66(3) 


Pure four dimensional site percolation |23 


0.161 


0.715 


0.683 


FKL97 pi 




0.50(4) 


0.61(4) 



The vector nature of the monopole current suggests that our model do not lie in the 
same universality class that the four-dimensional site percolation. 

We also see that our values are closed to the non-compact QED ]T2|], but more data is 
required to decide if the two models belong to the same universality class. 

Finally, we would like to stress the remarkable agreement with ||16|1 . Nevertheless, they 
perform their analysis using a different value of the Higgs coupling, (3 H ~ 0.9. As it pointed 
out in |[6|, a more careful analysis is necessary to decide if the critical exponents are the 
same along the percolation line beyond the end point. 



V. CONCLUSIONS 



We have performed a numerical simulation of the compact lattice gauge U(l) coupled to 
Higgs fields of unitary norm. Some results concerning the phase diagram previously obtained 
by Alonso et al. J0|, have been confirmed but measuring at the same time the behavior of 
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monopoles, condensation and percolation. The location of the percolation transition line 
-the line defined by the maxima of the monopole susceptibility- in the plane has 
been determined. This line overlaps the confined-Coulomb and the confined-Higgs phase 
transition lines, which are originated by a monopole-condensation mechanism, but it con- 
tinues away from the end-point. This is in contrast with the behavior of the monopole 
percolation transition in the non-compact abelian Higgs model where it is unrelated to any 
phase transition. In addition, the critical exponents of the monopole percolation transition 
in the region far away from the phase transition lines have been determined. The finite size 
scaling has been performed in terms of the monopole density instead of the coupling, because 
the density seems to be the natural parameter when dealing with percolation phenomena. 



While this paper was written up we have received Ref. |L6| where a model with both, 
scalar and fermion matter fields has been considered. Their results for the scalar sector are 
in perfect agreement with our simulations. It is interesting to compare the finite size scaling 
analysis for the critical exponents in the pure-percolation region obtained using couplings 
and density parameterizations. The values obtained in both cases are in perfect agreement. 
However our density parameterization results have been obtained with lattice sizes smaller 



than those in Ref. [16|. This suggests that the finite volume effects are smaller if our 



parameterization is used. 
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